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An elementary application of Algorithmic Complex- 
ity Theory to the polygonal approximations of curved 
billiards-integrable and chaotic-unveils the equivalence of 
this problem to the procedure of quantization of classical 
systems: the scaling relations for the average complexity of 
symbolic trajectories are formally the same as those gov- 
erning the semi-classical limit of quantum systems. Two 
cases-the circle, and the stadium-are examined in detail, 
and are presented as paradigms. 
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I. INTRODUCTION 

Chaos is certainly the most significant concept that 
has issued from the theory of dynamical systems and 
yet its true meaning, most concisely and universally 
encompassed in the equation chaos = deterministic 
randomness has not been fully adopted in the litera- 
ture and in the scientific community. This is some- 
how paradoxical, for even in popular magazines the 
idea has spread that chaos theory might be considered 
the third scientific revolution of the century-after rel- 
ativity and quantum mechanics. In my opinion, two 
are the main reasons behind this failure: firstly, the 
information-theoretical concepts implied in the notion 
of deterministic randomness are unfamiliar to most 
scientists; secondly, for the vast majority of physicists 
the "true" mechanics is not classical-where chaos is 
commonly found-but quantum, where chaos is, quite 
significantly as we shall see, absent. 

A clash is implied in this last statement: if chaos 
is absent in quantum mechanics, shouldn't it be also 
in classical mechanics, which is just the limiting case 
of the former ? This clash has lead people to draw 
all possible sorts of conclusions. Many have claimed 
that classical chaos must imply quantum chaos, via 
the correspondence principle [El, while others on the 
same basis have pretended that classical chaos theory 
should be derived from quantum dynamics. Others 



yet have argued that being chaos absent in the quan- 
tum mechanical theory of nature, it should be absent 
in nature altogether, and in particular at the macro- 
scopic level, and so farewell classical chaos. Taken 
at face value, this last statement implies a logical in- 
consistency, for quantum mechanics is just a theory, a 
description of nature. Yet, many physicists consider it 
a very good description in all respects, including chaos 
or its lack. The fact remains that chaos in nature is 
an undisputed reality, quite well described by classical 
dynamics. In previous works we have put forward the 
idea that the gap between the two mechanics, classical 
and quantum, is wider than can be naively expected 
from the correspondence principle |^-|| . 

In this paper I shall show that a similar situation is 
met in the dynamics of classical billiards: here, ratio- 
nal polygonal billiards play the role of quantum sys- 
tems, whose "classical limit" are curved billiards, to 
which they tend geometrically as the number of polyg- 
onal sides increases indefinitely. Indeed, while it has 
been since long recognized that the former are mod- 
els of non-chaotic behaviour [p|-H , one can use them 
to approximate chaotic curved billiard tables to ar- 
bitrary precision, and ask what happens then of the 
character of their motion. Even more: recently some 
observations of positive "effective" Lyapunov numbers 
in polygonal billiards have been published pd| , and a 
paradox of the same flavor as "quantum chaos" seems 
to arise. To resolve this paradox, in this paper I shall 
introduce an elementary, "physical" version of algo- 
rithmic complexity theory. In developing this theory, 
it will become clear that the procedure of approximat- 
ing curved billiards by polygons is quite analogous to 
that of quantizing classical systems: understanding 
the complexity of the motion in polygons can then 
be profitably used to clarify the issues involved in the 
other, more important problem. 

Our arguments are organized as follows: in the next 
section we review the fundamental pieces of informa- 
tion on integrable and chaotic billiards, and the no- 
tion of algorithmic integrability. To adapt algorithmic 
complexity to physical purposes, in Sect. Ill a simple 
coding of trajectories in billiards is introduced, which 
translates these into symbolic sequences. In Sect. IV 
the attention is focused on the scaling of the complex- 



ity of these sequences with respect to time, within 
certain time intervals: this leads to the concept of 
randomness (or order) within a range. The theory is 
immediately applied to the case at hand: Sect. V in- 
troduces two families of billiard tables: the circle (an 
integrable system), the stadium (a fully chaotic one), 
and their rational polygonal approximations, which, 
by analogy to quantum mechanics we also call polygo- 
nal quantizations. Rather than studying orbital com- 
plexity directly, I introduce here the concept of aver- 
age coding lengthy which has a clear physical meaning, 
and can be used to estimate the former. The case of 
circle quantizations is then studied in Sect. VI, where 
the paradox presented above is resolved. The stadium 
billiard is put to the same test in Sect. VII, and the 
problem of correspondence is addressed. Finally, I 
shall try to draw some general conclusions asking a 
challenging theoretical question. 

It is a pleasure and a honour for me to to dedicate 
these reflections to Boris Chirikov for his seventieth 
birthday as my teacher and colleague Joe Ford would 
have certainly done, had he been here now. 



II. CHAOS AND COMPLEXITY IN 
BILLIARDS 

Billiards are perhaps the dynamical systems that 
require less of an introduction, and just a few for- 
mal definitions arc necessary. A billiard table B is a, 
bounded, connected domain of the plane, with piece- 
wise smooth boundary. An ideal billiard in B is the 
dynamical system originating from the uniform mo- 
tion of a point particle-a ball-inside B, with elastic 
reflections at the boundary, following the familiar law 
of angle of incidence equal angle of reflection. 

In many games that can be played on a billiard, at- 
tention is paid at such reflections: one discretizes time 
(the integer n meaning the time of the n-th rebound) 
and considers the subset S of the tangent space of B, 
consisting of unit tangent vectors attached at bound- 
ary points, pointing inside B. S can be easily param- 
eterized by the pair (1,(1)), where I is the arc length 
along the boundary, and </> is the angle between the 
unit vector and the inner normal to the boundary, 
— f < (/) < f . In so doing, dynamics is a function T 
from S to itself, that maps the bounce occurring at 
In-i, with "exit angle" 4>n-i, into the new collision 
point In, and exit angle (/)„: 



{ln,4>n) = T{ln^i,(l)n-l). 



(1) 



Since this mapping preserves the canonical measure 
d^ = cos (j) dl d(j), billiards are among the simplest 



and most successful examples of conservative dynam- 
ics. But perhaps they mostly owe their success to the 
fact that a member of their family can be found at 
virtually all levels in the famous ergodic hierarchy: 
there are integrable ones-the circle-as well as K-ihe 
stadium, which is therefore also ergodic, and mixing. 
Donald Ornstein has shown (with Gallavotti [|2| ) that 
billiards can also be Bernoulli, and has conjectured 
(with Weiss [O]) that chaos in nature is mostly of 
this type: not without reason, we can say that bil- 
liards have served to shape our view of reality. In a 
sense, the present paper aims at the same ambitious 
goal. 

An interesting sub-class will be studied here, which 
covers part of the ergodic hierarchy, but falls short of 
producing chaotic representatives: the polygonal bil- 
liards 1^. Sinai has indeed proven g that these bil- 
liards have null metric entropy, and Ford has termed 
polygonal billiards with rational angles algorithmically 
integrable [A) W%, where the letter A was also in- 
tended to honor the memory of V.M. Alekseev, and 
his work [ll^ which showed that orbits of null entropy 
systems have null algorithmic complexity. The termi- 
nology makes it evident that a shift in perspective has 
taken place: while the ergodic hierarchy is concerned 
with statistical properties of ensembles of orbits, algo- 
rithmic theory deals with a new object, the complex- 
ity of the description of the motion, which will be the 
basis of our investigation p6[ . 

Seminal work |l^ on 74-integrable billiards is 
Zemlyakov and Katok's study of polygons whose 
vertex angles are all rational multiples of tt: it shows 
that these billiards satisfy the conditions for integra- 
bility except for the effect of vertices. Rationality of 
the angles provides a second constant of the motion 
(the angle of reflection </> times a suitable integer mul- 
tiple of 7r), but splitting of trajectories heading on a 
vertex provides the error |Tq] which prevents the sys- 
tem from being integrable. Notwithstanding these er- 
rors trajectories are still computable, in the sense that 
effective algorithms can be devised, in such a way that 
the number of informational bits in the output (the 
trajectory) is much greater than the number of infor- 
mational bits in the input (the algorithm plus the ini- 
tial condition of the motion) . Echardt et al in Q call 
such trajectories algorithmically meaningful. There- 
fore, ^-integrable systems are com put ationally akin to 
Liouville-Arnol'd integrable ones [|9[. We shall come 
back to these concepts in the following. 

A few years after ref. |^ Ford and other co-workers 
returned to the theme of rational billiards, presenting 
a seemingly different set of conclusions |11| : they ex- 
amined the rate of divergence of nearby trajectories 
and found that this rate is exponential, even for ra- 



tional billiards. The key factor in their derivation is 
the fact that nearby trajectories differ at time zero by 
a finite, fixed amount, and they are re-initialised to 
this fixed amount at each iteration of the Benettin- 
Strelcyn algorithm which computes Lyapunov num- 
bers. In the mind of these authors, this upper bound 
to precision stands for the human limitation against 
the dogma of infinite precision. If we take this limita- 
tion into account, they claim, trajectories effectively 
live on a multi-sheeted surface to which splitting at 
rational vertices gives an average negative curvature: 
the resulting motion is, practically speaking, chaotic. 
The first aim of this paper is to put order in these con- 
flicting observations, by utilizing a scaling approach to 
algorithmic complexity theory. But first, let us define 
the rules of our game. 



III. SYMBOLIC SEQUENCES IN BILLIARDS: 
A GAME 



The billiard in a circle (Fig. P is a noticeable exam- 
ple of Liouville-Arnol'd {L-A) integrability: the sec- 
ond, smooth integral of the motion being the angular 
momentum with respect to the center. Cutting the 
circle in two equal pieces, and inserting a rectangular 
strip in between the two halves gives rise to a fully 
chaotic billiard: the stadium [EG|. We shall get rid of 
all symmetries in this geometrical figure, and study 
the quarter stadium (Fig. 0). Let us now replace the 
circular sides in both billiards by a polygonal approxi- 
mation with equal sides: it is apparent that this can be 
done so to form a rational billiard. In both these cases, 
Vega, Uzer and Ford have found exponentially diver- 
gent trajectories, within their approximation scheme 
of course. In the following, I shall show that the two 
cases are different indeed. To do this, I first need to 
introduce a symbolic coding of this problem. 

How to code a dynamical object into a symbolic se- 
quence is something which follows either from physical 
insight, or practical convenience, or mathematical ef- 
ficacy. In our case, we shall elect to code trajectories 
according to the bounce coordinate I alone: for this, 
we shall assume that the circular parts of the bound- 
ary (or their polygonal approximations) are divided 
into a finite number, 5', of equal regions, each of which 
corresponds to a symbol, a, which can be taken to be 
a natural number from to S—1. We also decide that 
bounces on the straight segments of the stadium will 
not be registered. This coding is indicated in Figures 
|l| and g, for S* = 2. In a billiard game like this, we may 
think of putting detectors all around the boundaries, 
which are set off whenever the ball hits them. In the 
polygons, nearby sides can be connected to the same 



detector, so that the total number of output channels 
is S, also when M-the number of sides-is much larger 
than S, and is allowed to increase, while keeping S 
fixed, as we shall do in the following. The history of 
a trajectory is then coded in the record cri,cr2, . . • of 
boundary reflections. 

Now, let us play a funny game: the ball is initially 
set still at a fixed point, and a test player-chosen ap- 
propriately among our pals-can aim it by hitting it 
properly by the conventional cue. As the departing 
angle (f) is varied different trajectories are initiated, 
and different symbolic sequences {aj} are recorded. 
Let us now ask the player to do it twice, that is, to 
aim the ball a second time so that the sequence of 
bounces he had obtained in the first shot is repeated 
[ pl[ . It comes as no surprise that he will not be able 
to set the initial angle (j) to exactly the same value in 
both tries: the difference in this quantities causes the 
two resulting symbolic trajectories to agree only over 
a finite time-span. Trajectories of this kind are pic- 
tured in Figures || and |[ But how is this sport related 
to algorithmic complexity and to our problem ? 




FIG. 1. The initial portions of two trajectories in a cir- 
cular billiard table, characterized by slightly different ini- 
tial angles (fi. Their common initial position is at the bot- 
tom of the table. Also shown is a coding of rebounds with 
two symbols, a = (light portion of the boundary) and 
a = 1 (dark portion). 




FIG. 2. Same as Fig. 1, now for the quarter stadium 
billiard table. Notice that rebounds are coded only on the 
circular part of the boundary, by cr = (light thick portion 
of the boundary) and by a = 1 (dark thick portion). 



IV. DETERMINISTIC RANDOMNESS OF 
FINITE TRAJECTORIES 

In the light of the previous section, a bilhard can 
be thought of as a computer program, that is, a 
Turing machine capable of taking a binary sequence 
Pi, . . . ,Pi-the program, to which we shall momen- 
tarily return-as input, and producing a symbolic se- 
quence, (Ti, . . . , fJAT, the output. Recall now that the 
algorithmic complexity K{N) of a sequence {uj}, for 
j = 1, . . . , A'^ is defined very roughly as the length of 
the shortest computer program capable of outputting 
the sequence, and stopping afterwards |2^j2^. We can 
say that the complexity of a sequence is the length of 
its shortest definition: 



Pi, 



outputs 

,PL => CTi, 



, (Tat implies K{N) < L. 



(2) 



Algorithmic complexity theory teaches us that most 
sequences-in probabilistic sense-of length N have 
complexity close to maximal, i.e. N. Moreover, in 
the limit of increasing N , almost all of them are ran- 
dom^ in the sense that their complexity K{N) grows 
as N |24|] . At the same time, computable sequences 
exist, for which K{N) is much less than iV, and grows 
less than linearly. We apply now this theory to the 



dynamical sequences constructed in the previous sec- 
tion. Let us call their complexity K(N^(j)), where we 
have explicitly indicated that this number may de- 
pend on the initial condition of the motion. We shall 
then follow Aleksecv (l|], Chirikov |2|l, and Ford (26|, 
and identify order with computability, and chaos with 
randomness. 

The central issue is then to find a computer code 
to output the desired sequence. For our dynamical 
sequences, cti, . . . , crjv, this can be obtained in any ab- 
stract language by an encoding of: a) the geometrical 
rules of the game, which require a fixed number of 
bits, Cmachino (which depends only on the machine on 
which the rules are coded), plus b) the instructions set 
for fixing the billiard boundaries, of length Cboundaryj 
plus c) the specification of N and d) of a certain num- 
ber of digits of (j). Accordingly, the length L{N, (/>) of 
this program can be estimated as 



L{N, (j)) ~ Cmachino + Cboundary + loga N + A{N, ( 



(3) 



where the function A(iV, 0) is defined as the number 
of bits of necessary and sufficient to determine the 
first N symbols in the sequence a. This function can 
serve to estimate the complexity K{N, (p) in its most 
important aspect: the N dependence. For it is clear 
that the opposed behaviors presented above-random 
and computable-are of deep physical significance: if 
in a range Ni < N < N^ the function A(A^, (j)) grows 
linearly, A(A^, (f) ~ \N, the complexity K{N, (f) has 
the same leading behaviour, and the sequence {ctj} 
should duly be termed random, or chaotic. When, 
in a similar interval, A(A^, (j)) grows less than linearly, 
{dj} should be called computable, or ordered. We 
therefore introduce the notion of order (and random- 
ness) within a range. The physical significance of such 
order (or randomness) will then be proportional to the 
importance of such range. 

We have emphasized this scaling approach in B to 
rebut a common objection to the application of com- 
plexity theory to finite dynamical sequences [M . The 
rationale behind this idea is evident, and is quite simi- 
lar to that adopted in the independent theory of com- 
putational time complexity: not so much the value of 
complexity is relevant, but the way it increases as the 
problem size grows, for this may render it quickly un- 
feasible. In fact, there is no point in practicing for our 
billiard player if he is playing a chaotic billiard: for 
each additional bounce he wants to set correctly his 
aiming precision in the initial angle (j) must increase 
geometrically. Out of metaphor, when the billiard is 
a system whose symbolic dynamics we want to pre- 
dict, to obtain a linear increase in forecast precision 



we must exponentially increase the accuracy in the ini- 
tial conditions [ psj . In most instances this demands 
an exponential increase of economical resources. 



V. THE CIRCLE AND THE 
STADIUM-POLYGONAL QUANTIZATION 

To draw a parallel with quantum mechanics |29|| 
suppose now that Nature, by means of our observa- 
tions, clearly shows that curved billiard boundaries 
are a mathematical idealization and that physically 
we can just have polygonal billiards with an arbitrary 
but finite number of sides, M . These sides will be 
inscribed in the circle and in the stadium, and I shall 
call the resulting polygonal billiards a quantization of 
the curved ones. In this context, M is a crucial quan- 
tity, which plays the role of a semi-classical parameter. 
Clearly, when we let M go to infinity we regain geo- 
metrically the original table. The question is, will we 
obtain the same dynamics ? 

To answer this question in a complete, algorithmic 
sense, let us consider again the coding introduced in 
Sect. Ill, and let us study the function A(A^, 0) de- 
fined in Sect. IV, or rather its inverse, Af(A, 0), which 
quantifies the number N of symbols Uj which can be 
predicted knowing A digits of 0. It is clear that order 
(or randomness) of a dynamical sequence can also be 
inferred from the study of N{A,(I)). Moreover, this 
function has a transparent physical meaning: if we let 
e :— 2~^, in this new independent variable the func- 
tion N{e, 4>) := N{— log2 £, (f>) can be called the length 
of the codable trajectory under uncertainty e. Let us 
go back to the player example. He is only capable of 
setting the angle with a certain error e. Under these 
circumstances, only finitely many dynamical symbols 
will typically coincide in the initial portions of a series 
of two shots a and a': their number is A^(A, (p), 

N{A,(/)) = k if ai^a[,...,(Jk^ <^k'<^k+i 7^ (r'k+i- 



(4) 



Following |30) we also say that N{e, (/)) is the first error 
time, when any two initially e-close trajectories start 
to have different symbolic sequences. 

Numerically this function is rather straightforward 
to compute, and its theoretical analysis can be carried 
out in full detail (see the Appendix and below). It is 
plotted in Fig. M versus 4> at fixed uncertainty e, for 
the M = 8 quantization of the circle (a regular octag- 
onal billiard). One observes that low values, like that 
seen at (j>i, are associated with trajectories heading 
very early in history towards a vertex, where a sym- 
bolic error may occur |pQ;p[. Large values of -/V(A, (j>) 



are found when this happens much later: 02 is associ- 
ated with a periodic trajectory which stays far away 
from the vertices ]3l| ]. 

It can be shown that the behaviour of N at the pe- 
riodicity values is at the root of the dynamical prop- 
erties of the system. In fact, let the ideal, frictionless 
ball run for an infinite time. A coding function F of 
the initial angle (p can be defined as 



n</.) = 5^a,5- 



(5) 



This function represents the translation dictionary be- 
tween the trivial code-the infinite sequence of digits 
of (f), and the dynamical code-the infinite sequences 
of bounces. In chaotic systems, like the motion 
over surfaces of constant negative curvature, and the 
anisotropic Kepler problem, the relation between the 
properties of N and the function F can be fully ex- 
posed. This relation justifies the multi-fractal prop- 
erties of the coding function F, which were originally 
investigated by Gutzwiller and Mandelbrot [p2|,P3| 
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FIG. 3. Coding length A''(A, 0) vs cp for an octagonal 
billiard, with 5* = 8. The starting point is at the center 
of one of the sides. Two values of e are used: e = .01 
(bottom continuous line) and e — .005 (top broken line). 
The arrows mark the values 0i and 02 described in the 
text. 

In player's terms, the angle 02 in Fig. @ is an easy 
shot, and 0i a tough one: we find then convenient 
to average over initial angles, so defining the average 
coding length A{e): 



Aie):^ J Ni^log^e, 



C0S( 



(6) 



which will become our main indicator of the complex- 
ity of the motion. A canonical average over the full 
phase space can be also defined and computed, with 



quantitatively similar results to those we are going to 
expose. Let us pause no more, and compare the be- 
havior of this quantity in the circle, the stadium, and 
their polygonal quantizations. 



VI. ORBITAL COMPLEXITY IN 
INTEGRABLE AND POLYGONAL BILLIARDS 

We have now set the stage for the examination of 
the average orbital complexity in our billiards. We 
start from the case of the circle and its polygonal ap- 
proximations: before resorting to exact analysis, let us 
have a look at the numerical data. Fig. plots A{s) 
versus e. The scale is doubly logarithmic: according 
to eq. (^ the power-law behavior of these curves for 
small £ is a manifestation of the ordered character of 
the motion, which is well assessed in the literature |g] , 

ra, a. 
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FIG. 4. Average coding length Am{s) vs e for the circle 
and its polygonal quantizations, with S — 2. The values of 
M are chosen to be twice a prime, and roughly geometri- 
cally increasing. The continuous line is the second formula 
in eq. (0) with M = 62. 

Yet, the physical picture evident in Fig. is much 
richer. Firstly, as M grows at fixed (large) e, AmIs) 
tends to A^c (e) , where the label M is self-explanatory, 
the earlier the larger the value of e. This is indeed 
clear: at large e the coding length is short, and differ- 
ences between polygonal boundaries with large M and 
the circle are not significant. But if we turn now our 
attention to the left part of Fig. H, where relatively 
smaller values of e are plotted, a less expected phe- 
nomenon appears: increasing M at fixed e the differ- 
ence between Am (s) and A^o (e) seems to grow rather 
than vanish. Moreover, this discrepancy is brought 
about by a decreasing Am{£)- Recall that Am{£) is 



the average of the inverse function of A(A^, (j)) (if you 
rotate the graph clockwise by ninety degrees, N ap- 
pears plotted on the horizontal axis and e — 2^^*^^^ 
on the vertical): in this region the complexity of a 
trajectory of fixed length N — Am grows when M 
is increased. This is the phenomenon observed by 
Vega, Uzer and Ford: polygonal billiards have their 
own brand of instability, generated by splitting of tra- 
jectories at vertices. The more the vertices, the larger 
the Lyapunov numbers computed in [ pT[ . The para- 
dox presented in Sect. II is so exposed. 

But we are now equipped to resolve it: observe first 
that the increase in complexity obtained by raising 
M at fixed e is only temporary: if we keep going, we 
find that Am{e) finds a minimum, and then inverts its 
course to reach quickly (we shall come back on the rate 
of this convergence later on) the circle value Aqo |34| . 
Furthermore, a deeper argument must be made: sit- 
ting at fixed e is not a proper thing to do, not even in 
the presence of human limitations to finite precision. 
In fact, we cannot increase precision indefinitely, but 
we certainly can over a finite, physically reasonable 
range. For instance, this may be dictated by the com- 
putational power of the machine on which Fig. has 
been computed. Exploring this range, from larger to 
smaller uncertainties, we discover that Am{£) starts 
off like the L-^-integrable circle, Aoo(e) ^ s^2 ^ and 
successively redirects its course on a different line, 
with the same exponent: a transition from L-A to 
^-integrability has taken place. 

These observations can be derived rigorously from 
the analysis developed in the Appendix, with the re- 
sult: 



AM{e) 



1 

/2e 



/Me 



for Me > H, 
for Me-^H, 



(7) 



where the crucial functional dependence on M is ap- 
parent, and where iJ is a constant, which plays a sim- 
ilar role to Planck's constant in the "usual" quantum 
mechanics. In the above equation, the first behaviour 
is that of the pure circle (correspondence). The sec- 
ond equation in (^ leads to an estimate for the av- 
erage complexity which increases only logarithmically 
with N and M: 



A{N) <C + 2 log2 N + log2 Af , 



(8) 



where C is a positive constant, independent of M and 
TV, and therefore from eq. (^ one obtains [Q 



L{N) <C' + 3 log2 N + 2 log2 M, 



(9) 



where C" is another positive constant. This result 
should be compared with the existing literature on 
the topological complexity of symbolic dynamics p& . 



In conclusion, we see that vertices add logarithmi- 
cally to the complexity, but only as long as Me ^ 7J, 
i.e. log A/ < logiJ + A. We have so reconciled the 
observation of rcf. ||ll| and the established knowledge 
on polygonal billiards. We can now turn to a more 
delicate problem: what happens of orbital complexity 
in the polygonal approximations of a chaotic billiard. 



VII. ORBITAL COMPLEXITY AND 

CORRESPONDENCE IN CHAOTIC 

BILLIARDS 

In the previous section we have established that ra- 
tional approximations of the circle cannot be called 
chaotic, not even under the assumption of finite pre- 
cision. Is this achievement possible to their more so- 
phisticated relatives, those inscribed in the stadium ? 
After all, as M tends to infinity, these fellows tend 
to a fully chaotic system, and one thing is certain, 
correspondence principle must be obeyed |39[|. 
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FIG. 5. Average coding length Am{£) vs e for the quar- 
ter stadium and its polygonal quantizations, with S = 2. 
The values of M are powers of two. The curves l(x) and 
h{x) are given in eq. (FlQ) with M = oo and M = 128, 
respectively, and A = -1.731, B = 0.831, p = 1.78. 

Fig. H is the analogue of Fig. H, now for the quarter 
stadium billiard table. The abscissa is e in log scale, 
i.e. minus A. The logarithmic character of the curve 
AoQ clearly reveals that the complexity of trajecto- 
ries grows linearly with length: chaos is here mani- 
fest in its essence. For large values of Ale, Am{£) is 
roughly equal to Aoo(e). For instance, over the inter- 
val (10~^, 10~^) ^32768 is a logarithmic function of e; 
we can therefore expect that a corresponding range in 
N exist so that the average complexity K{N) grows 
linearly: here, trajectories are random, in the sense 



explained in Sect. IV. Notice that the existence of 
this range is what permits "sensible" numerical ex- 
periments of chaotic motion. 

The effect of vertices noticed in the circle quanti- 
zation appears again when following the curve to the 
left, as it leaves Aoo staying slightly lower than this 
latter, i.e. showing a relative increase in the com- 
plexity of the motion. As noted above, this increase 
is quantified by the logarithmic contribution logM. 
Alas, this excess of zeal rapidly turns into failure: 
the algorithmic simplicity of polygonal boundaries is 
soon detected, and A32768 starts to grow like e~2. 
What has happened is that the complexity estimate 
(H), previously dominated by that of the curved sta- 
dium, takes over for good. In fact, an exact analysis 
can be performed here too, showing that 
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Me > H, 
Me<^H, 



(10) 



where A, B, p and H are suitable constants, the last 
playing the same role as in the previous section. Con- 
sequently to this equation, in the region Ale <ti H 
complexity estimates of the form (g), (g) still apply. 

The end of this paper is in sight, and we must start 
drawing conclusions. There is a significant difference 
between the case of the circle (Fig. ^ and eq. (0) 
) and of the stadium quantization (Fig. ^ and eq. 
([lO|)): while the former, as we noted, shows a transi- 
tion from A to L-^-integrability, the second matches 
^-integrability with randomness, and this mating is 
troublesome, to say the least. Perhaps the most evi- 
dent consequence is the following: let me ask how long 
must a symbolic sequence be to perceive the finitc- 
ness of the number of boundary sides via its algo- 
rithmic manifestations. The answer is instructive. In 
the polygonal quantization of the circle, and of the 
stadium as well, the critical value eM at which the 
AI = 00 curve is abandoned scales as £m ~ ^- Yet, 
at the same time the different structures of the two 
problems imply that the length of polygonal-curved 
accordance is proportional to vM in the integrable 
case, and only to the logarithm of M for the chaotic 
stadium ||§. 

These results have been obtained by explicit calcu- 
lation. Yet, I believe their algorithmic nature to be 
more important: we must realize that both circle and 
stadium polygonal quantizations are dynamical sys- 
tems endowed with an amount of complexity which 
scales as the logarithm of M. Yet, while the M-circle 
dynamics unfolds this conrplexity slowly, at a loga- 
rithmic pace, its M-stadium relative does it eagerly, 
linearly in time, so that the complexity reservoir is 
quickly exhausted. In the Conclusions I shall briefly 



present my views on the physical imphcations of these 
facts. 



VIII. CONCLUSIONS: CHAOS IN NATURE, IS 
THERE ANY ? 



Algorithmic complexity theory, used here in a physi- 
cist's fashion via the concept of average first error 
time, clarifies the simplicity of the dynamics of polygo- 
nal billiards, and offers a solution of the paradox which 
arises from the juxtaposition of ref. M and O] . At the 
same time, this theory provides a complete character- 
ization, in all ranges of parameters, of the dynamical 
properties of the families of billiards studied in this 
paper. These latter, in turn, are significative repre- 
sentatives of integrable and chaotic billiards, and we 
can safely expect our results to be quite general. 

In essence, the algorithmic estimates derived in this 
paper are manifestations of the same phenomenon ob- 
served in the Schrodinger quantization of the Arnol'd 
cat i|,|,|,| (see eqs. (8) and (9) of ref. ||), of 
bounded systems B, and in the classical dynamics 
of discrete systems (see Sect IV of ref. Q). Truly 
then, the nature of chaos and order is information con- 
tent ! Yet, if we admit this, we must be prepared to 
bear the consequences when going back to the problem 
of quantum mechanics: then, the metaphor of ratio- 
nal billiards exposed here permits to understand the 
claim that the correspondence principle is validated 
for integrable systems, and violated for chaotic ones. 
It was shown in the last section, in fact, that in the 
first case a discrete nature-in which only polygons are 
allowed-would let us play with our curved theoretical 
models for a long time. For the same reason, in the 
quantum mechanics of an integrable system, the ac- 
tion of increasing a semi-classical parameter provides 
a complexity reservoir which keeps up with the classi- 
cal description for a time-span which grows apprecia- 
bly, as a power-law, in this parameter. In the second 
case, to the contrary, the time of chaotic freedom is 
logarithmically short, and the essence of correspon- 
dence, to regain classical/curved complexity by quan- 
tum/polygonal computations is exponentially remote 
in the semi-classical parameter pl| . Quantum dynam- 
ics is a computable theory, and can "correspond" to 
the uncomputable classical mechanics only as far as 
its algorithmically simple nature allows it |4^ . 

Notwithstanding this evidence, people have been for 
a long time reluctant to admit that there is not such a 
thing as "quantum chaos" , or that this chaos is some- 
thing different-and less-than deterministic random- 
ness. This hindrance is now history. To the contrary, 
Boris Chirikov has appreciated this fact since the very 



beginning, and his notions of transient chaos |g^ and 
pseudo-chaos Q , [Q , show this very clearly. More- 
over, he has put a strong emphasis on this latter con- 
cept: largely simplifying, but I believe appropriately, 
one could say that Chirikov views the pseudo-chaos 
(that we call ^-integrability) of discrete systems (and 
by consequence of the digital computers as well) as a 
faithful representation of what takes place in quantum 
mechanics pq]-and in nature, at all levels. At this 
very last point we part company: in fact, an unset- 
tled debate still holds on the meaning and validity of 
the correspondence principle p7[ , and on the question 
whether quantum ^-integrability is enough to cope 
with a world in which chaos seems to be essential. 
In my works with Joe Ford a negative answer to this 
question has been presented, suggesting that-in this 
sense-quantum mechanics might be incomplete^ and 
that the "reductionist" program of deriving classical 
chaos from quantum mechanics could not be achieved 

Certainly, classical mechanics is inadequate to de- 
scribe reality, but it contains the gene of non- 
computability, which we believe should be present in a 
complete theory of nature; at the same time, the more 
fundamental quantum mechanics inherited this gene 
in too a tame form to be effective. This is, in nuce, 
the point of contention. Will time settle the matter- 
or will it be that the general indeterminacy principle 
foreseen by Ford p9] will change everything around ? 



IX. APPENDIX: AVERAGE FIRST ERROR 
TIME CALCULATION 

We shall now justify the formulae presented in 
Sects. VI and VII via a direct calculation of the aver- 
age first error time A(e). Firstly, let us consider the 
case of the (integrable) circular billiard. The angle (p 
is a constant of the motion: letting u := t: — 2(p the 
motion is Z„ — /„_i +uj, the usual ergodic (for irra- 
tional w/tt) rotation of the circle. Trajectories with 
different initial angle separate at a linear speed 2e, 
where e has the same meaning as in the main body 
of the paper (see Fig. |l|). Consequently, a bunch 
of trajectories opens up linearly in time like a slightly 
unfocused light beam. An error occurs when a bound- 
ary point of the circle partition which determines the 
coding enters this light cone. Let us now estimate the 
average time required for this to happen. 

Let S to be the number of cells of the partition of 
the boundary, and Ss '■— -S- their common length. Let 
also 0n '■— In rnod Ss, n = 1,2, . . .. According what 
we have just said, an error occurs when the angle 0„ 
of the reference trajectory (the center of the beam) 



comes within ne/2 of zero or one, where e := j-e. 
Let Pn(e) be the probabihty that the first error occurs 
at time n. Ergodicity impHes that pi{s) = e. Certain 
approximations are required to evaluate p„ for n > 1: 
we assume that {On}, n = 1,2,... are uncorrelated 
random variables, so that p„(e) — "-ell^Ji (1 ~ J^)- 
Estimating then the asymptotic (small e) behaviour 
of A{e) = J2n^Pn{£) by standard techniques gives 
A{£) ~ -v/f^"^' ^^^ therefore 



A{s) 



-s- 



(11) 



which well reproduces the behaviour numerically 
found in Sect. VI. 

Surprisingly, the same e dependence can be shown 
to hold for the polygonal billiards considered in this 
paper. The presence of the vertices imposes a more 
sophisticated analysis. When the beam bounces fully 
on the same polygonal side, its angular amplitude is 
conserved, and its transverse dimension grows linearly 
in time on average, exactly as it does in the circular 
billiard. To the contrary, when the beam impinges 
on a vertex of amplitude a, its angular amplitude is 
increased by the amount 2(7r — a). In the polygonal 
approximations of the circle, we have a — {1 — 2/M)tt, 
while a = (1 — l/2A/)7r for the quarter stadium, so 
that (/)^j := 4t is the resulting perturbation in the 

(s) 

circular case, and (j)\J :— p- in the quarter stadium. 

Two regimes must now be considered. Firstly, when 
e is much larger than (f>M, the effect of vertices is 
negligible, trajectories behave as in the correspond- 
ing curved billiards, and AM(e) coincides to a good 
approximation with Aoo(e). In the opposite regime, 
when (f>M is much larger than e, impinging on a vertex 
causes a major enlargement of the beam. We can as- 
sume that this leads almost immediately to an error: 
a similar analysis to that performed in the circular 
case can then be carried out. We let 6m = ^ be the 
approximate length of the polygonal sides: this quan- 
tity plays here the same role 6s did above. The result 
is then 
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for the Af -circle, and 
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(12) 



(13) 



for the Af -stadium, where p is a parameter which mea- 
sures the effective length of trajectories in between two 
coded bounces. Finally, the intermediate region be- 
tween the two regimes is determined by the condition 
(j)M ~ £ described above, which becomes the "indeter- 
minacy principle" Me ^ H , where H is a constant. 
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